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ABSTRACT 

We have carried out high resolution observations with Atacama Large Millime- 
ter/Submillimeter Array (ALMA) of continuum emission from Orion KL region. We 
identify 11 compact sources at ALMA band 6 (245 GHz) and band 7 (339 GHz), in¬ 
cluding Hot Core, Compact Ridge, SMA1, IRc4, IRc7, and a radio source I (Source I). 
Spectral energy distribution (SED) of each source is determined by using previous 3 mm 
continuum emission data. Physical properties such as size, mass, hydrogen number 
density and column density are discussed based on the dust graybody SED. Among 
11 identified sources, Source I, a massive protostar candidate, is a dominant energy 
source in Orion KL. We extensively investigate its SED from centimeter to submillime- 
ter wavelengths. The SED of Source I can be fitted with a single power-law index of 
1.97 suggesting an optically thick emission. We employ the H~ free-free emission as 
an opacity source of this optically thick emission. The temperature, density, and mass 
of the circumstellar disk associated with Source I are constrained by the SED of H~ 
free-free emission. Still the fitting result shows a significant deviation from the observed 
flux densities. Combined with the thermal dust graybody SED to explain excess emis¬ 
sion at higher frequency, a smaller power-law index of 1.60 for the H~ free-free emission 
is obtained in the SED fitting. The power-law index smaller than 2 would suggest a 
compact source size or a clumpy structure unresolved with the present study. Future 
higher resolution observations with ALMA are essential to reveal more detailed spatial 
structure and physical properties of Source I. 

Subject headings: ISM: individual objects (Orion KL) — radio continuum: stars - 
stars: formation — stars: individual (Orion Source I) 
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Introduction 


Located at a distance of 420 pc fro m the Sun (IHirota et alJl2007l : lMenten et al.ll2007l : lKim et al 


2008 ), the Orion Kleinmann-Low (KL; K leinmann &; Low 1967 ) region is known as the nearest site 
of active massive-star formation. It has been recognized as one of th e best laboratories to study 
massive star-formation processes (IGenzel fc Stutzkil Il989l : iBallvl 120081 ) . Because of its extremely 
high opacity, observations have been made mainly in the centimeter, millimeter/submillimeter, and 
infrared wavelengths. In particular, radio interferometers have been powerful tools to s tudy in 


detai l about physical properties of young stellar objects (YSOs) in Orion KL (e.g. iBeuther 


Plambeck et al 


2004, 

20(3, 

2006; 

Tang et al. 

2010; 

Favre et al. 

2011; 

Friedel & Weaver 

2011; 

Zapata et al. 

2011; 


eta 


2013l l at the highest spatial resolution. These observ ational studies have identified 


some of the remarkable sources, such as Becklin-Neugebauer (BN; Becklin & Neuge baue r 1967) 
object, Hot Core, Compact Ridge, a radio source labeled I (Source I), infrared source labeled n 
(Source n), and submillimeter source identified by Submillimeter Array (SMA), SMAl, although 
their properties are still a matter of debate. 

Among th ese co mpa ct sou rces, the most prominent energy source in Orion KL is thought to be 


Source I ( Ment en fc Reid 19951 ). It is a candidate of massive protostar driving a so-called low veloc- 


thermal SiO lines and H?0 masers 

Genzel & Stutzki 

1989; 

Gaume et al. 

1998 

; f 

lirota et al. 

2007; 

Plambeck et al. 2009; ZaData et al. 

2012; 

Niederhofer et al. 

2012; 

Greenhill et al 

2013 

; Kim et al. 

2014 

). At the center of this outflow, there is a cluster of vibrationally excited SiO masers tracing 


a disk wind emanating from the surface of the circumstella r disk with a diameter of ~100 AU 


(Kim et al. 


2008 


Matthews et al.ll201(l : ICreenhill et al .11201.4 1. A compact radio continuum source 


is as sociated at the center of the SiO masers which is inter preted as an edge-on circumstellar 


disk (IReid et al.ll2007l : ICoddi et al 


2011 


Plambeck et alJl201.ll ). Such a disk-jet system is recently 


studied in the submillimeter H 2 O line s by using a newly constructed Atacama Large Millime¬ 
ter/Submillimeter Array (ALMA) (Hirota et ah 2012, 2014a). 


Nevertheless, the nature of Source I and associated disk/outflow system are still far from a 
complete understanding. One of the long-standing issues is an origin of the radio continuum emis¬ 
sion associated with Source I. A spectral energy distribution (SED) of Source I from centimeter 
to submillimeter wavelength can be fitted to a power-law function, F„ oc v 2 , indicative of opti¬ 


cally thick emission up to submillimeter wavelength (Plambeck et al. 2013, and references therein). 
Based on such a SED of Source I, it has been interpreted as either a combination of free-free emis - 


2006 


Reid ct al 


2007; 


sion and thermal dust emission (IBeuther et al- 2004, 

or a free-free emission from a neutral atomic/molecular gas (lReid et al 


Plambeck et al 


2007; 


2013) 


Plambeck et al 


2013). Although physical properties of Source I such as mass, temperature, and ionization degree 
are still under debate, more recent results support a latter scenario (e.g. P lambeck et al. 201 3j). 


In order to reveal physical properties of continuum sources in Orion KL, in particular Source I, 
we present observational study of multi-band continuum emission with the sub-arcsecond resolution 
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by using ALMA. The ALMA observations open new wavelength windows with the higher spatial 
resolution. Combined with previous high-resolution observations at lower frequency bands, we will 
discuss about basic physical properties of some of the key sources in Orion KL. 


2. Observations and Data Analysis 

Observations of the millimeter and submillimeter continuum emission were carried out with 
ALMA in several sessions during the early science operation in the cycle 0 period at bands 6 and 7 
(ADS/JAO.ALMA#2011.0.00199.S). We also employed the ALMA Science Verification (SV) data 
(ADS/JAO.ALMA#2011.0.00009.SV) at band 6. Details of each session is summarized in Tablc[l| 


2.1. Band 6 data 


The observation in cycle 0 was done in the extended configuration on 2012 April 08 with 
17x12 m antennas. The primary beam size of each 12 m antenna is 25 ,/ . The baseline lengths 
ranged from 17 to 310 kA (from 21 to 385 m). The observed frequency ranges were 240—244 GHz 
and 256—260 GHz consisted of 4 spectral windows with bandwidth of 2 GHz for each. Dual po¬ 
larization data were obtained simultaneously. The ALMA correlator was set for low resolution 
wideband continuum observations and the spectral resolution was 15.625 MHz. The target source 
was Orion KL and the tracking center position was taken t o be t he bur sting 22 GHz H 2 O maser, 


RA(J2000)=05h35ml4.1250s, Decl(J2000)=-05d22'36.486"(Hirota e_t al. 201l|, 2014a). The total 
on-source integration time was 30 s. A primary flux calibrator, band-pass calibrator, and sec¬ 
ondary gain calibrator were Callisto, J053851-440507, and J0607-085, respectively. The system 
noise temperatures ranged from 102 K to 122 K depending on the spectral window. 


To improve the image quality, we combined the ALMA SV data at band 6 with our cycle 0 
data. The SV data were obtained in the 20 frequency settings to cover full spectral emissions from 
214 GHz to 247 GHz. We here only employed part of the frequency range from 230 to 232 GHz 
to achieve an image sensitivity comparable to the cycle 0 data. The SV data were taken on 2012 
January 20 in the compact configuration with 16x12 m antennas separated from 14 to 203 kA (from 
17 to 265 m). The tracking center position of Orion KL was set to be R. A. =05h35ml4.35s and 
decl.=—05°22'35".0 (J2000), which is 4" northeast of that of the cycle 0 data. 

In this paper, we define the band 6 frequency to be 245 GHz as it is the central frequency of 
all the combined data. 
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2.2. Band 7 data 


For band 7, observations in two different frequency settings were carried out; continuum 
mode and spectral line mode. The tracking center position was the same as band 6 cycle 0 data; 
RA(J2000)=05h35ml4.1250s, Decl(J2000)=-05d22'36.486". 

Observations in the continuum mode was done in the extended configuration on 2012 Oc¬ 
tober 23 with 23x12 m antennas. The primary beam size of each 12 m antenna is 17". The 
baseline lengths ranged from 19 to 425 kA (from 17 to 372 nr). The observed frequency ranges 
were 341.5—345.4 GHz and 353.5—357.5 GHz consisted of four 2 GHz-bandwidth spectral windows 
with dual-polarization. The ALMA correlator was set for low resolution wideband continuum ob¬ 
servations and the spectral resolution was 15.625 MHz. The total on-source integration time was 
100 s. A primary flux calibrator, band-pass calibrator, and secondary gain calibrator were Callisto, 
J0423-013, and J0607-085, respectively. The system noise temperatures ranged from 116 K to 
162 K depending on the spectral window. 


We carried out monitoring observations of the submillimeter H 2 O lines at band 7 to study the 
H 2 O maser burst event at 22 GHz in multi-frequency ( Hiro ta et a l. 2014a ). For this purpose, three 
sessions of spectral line observations were carried out on 2012 July 16, August 25, and October 21 
during the cycle 0 period. Number of 12 m antennas were different from epoch to epoch; 21, 28, and 
22 at the first, second, and third epochs, respectively. Array configurations were in the extended 
mode with the baseline lengths of 16-427 kA (15-398 nr), 22-416 kA (21-388 nr), and 18-389 kA 
(17-363 nr) at the first, second, and third epoch, respectively. Four spectral windows were set at 
the frequency of 321.0-321.5 GHz, 322.1-322.6 GHz, 334.4-334.9 GHz, and 336.0-336.5 GHz with 
the bandwidth of 469 MHz for each. Spectral resolution was set to be 0.122 MHz. The on-source 
integration time of the target source was about 100 s for each session. A primary flux calibrator, 
band-pass calibrator, and secondary gain calibrator were Callisto, J053851-440507/J0423-013, and 
J0607-085, respectively. The system noise temperature ranged from 125-341 K, at which frequency 
around 322 GHz was significantly affected by the strong atmospheric absorption of the 325 GHz 
H 2 O line. 


Similar to the band 6 data, we define the band 7 frequency to be 339 GHz as the central 
frequency of all the combined data. 


2.3. Synthesis imaging 

The data were calibrated and imaged with the Common Astronomy Software Applications 
(CASA) package. For the band 6 data, we combined calibrated SV and cycle 0 data by using 
the task cone at in CASA. Plot of uv coverages of full-array are shown in Figure [1] For band 7 
data, we combined all of the calibrated data for both in continuum and spectral line mode by 
using the CASA task concat. After combining the data, continuum images of Orion KL were 





Table 1. Summary of Observations 


Band 

Date (in 2012) 

Center 

frequencies 51 

[GHz] 

Total 

bandwidth 

[MHz] 

Effective 

bandwidth 

[MHz] 

Number 

of 

Antennas 

On-source 

Time 

[sec] 

FWHM 

[arcsec] 

PA 

[degree] 

rms 

[mjy beam - 

- 1 ] 

6 

Jan. 20 

231 

1875 

167 

16 

1196.1 

1.76x1.19 

-i 


19 


Apr. 08 

250 

8000 

3547 

17 

30.3 

0.68x0.50 

-75 


5 


All uv 

245 





0.81x0.64 

-70 


9 


uv>100 kA 

245 





0.63x0.47 

-72 


5 

7 

Oct. 23 

350 

8000 

3641 

23 

100.8 

0.47x0.44 

0 


9 


Jul. 16, Aug. 25, Oct. 21 b 

329 

1875 

648 

21,28,22 

326.8 

0.55x0.46 

59 


12 


All uv 

339 





0.49x0.45 

59 


9 


uv>100 kA 

339 





0.46x0.41 

57 


5 


a Center frequencies for LSB (lower side band) and USB (upper side band). 

Only one spectral window in LSB is employed for band 6 data on January 20, 2012 (SV data). 

b Combine three epochs for spectral line monitoring at band 7. 
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made using the CASA task clean. The line emissions were excluded to make continuum images 
by integrating over the line-free channels. Resultant effective bandwidths were almost half of the 
observed frequency range as listed in Table [I] Both phase and amplitude self-calibrations were 
done with the continuum images by the CASA task gaincal and these results were applied to the 
visibility data by the CASA task applycal. The resultant image mis noise level of the continuum 
emission and the uniform-weighted synthesized beam size were summarized in Table |l] 


3. Results 


Figures [2] and [3] show continuum images of Orion KL at bands 6 and 7, respectively, obtained 
with various uv coverages. To achieve higher image quality and sensitivity, we combined different 
datasets observed in different uv coverage as explained in previous section. Overall spatial structures 
of the band 6 and band 7 images with full uv coverage (Figures Oja) and0Ka)) appear quite similar 
with each othei\_ The main emission region shows an elongated ridge-like structure (“main dust 
ridge” in Tan g et al. 2010j) along the northeast-southwest direction consisted of Source I, Hot Core 


and SMA1 as labeled in Figure 01 At western side of this main dust ridge and south of the BN 
object, there are several diffuse emission components located along the north-south direction, which 
can be identified as the Northwest Clump and Compact Ridge (Figure 01). Our ALMA images trace 
part of these extended emission corresponding to compact highest density condensations. 


In order to identify compact continuum sources, we made synthesized images by employing 
only baselines longer than 100 kA as shown in Figures 0Kb) and 0Kb). The extended emission 
components are resolved out in our ALMA images (Figures 0Kb) and 0Kb)) because of the lack 
of short baselines in the ALMA extended configuration. Although the high resolutions images 
(Figures 0Jb) and 0Kb)) show lower intensities than those of full uv coverages (Figures 0K a ) an d 
0Ka)), compact sources are clearly detected with sufficiently high signal-to-noise ratios greater than 
5a. Note that there are negative contour levels in these images due to insufficient uv coverages 
for shorter baselines. When we employ the visibility data with baseli nes longer than 200 kA, only 
Source I remains unresolved implying a compact structure dBeuther et al. 2004). 


For band 6 data, we compare the cycle 0 data (Figure0Kc)) and SV data (Figure0Kd)) obtained 
in the extended and compact configurations, respectively. Only the higher resolution cycle 0 data 
resolve the compact cores, as can be seen in the main dust ridge (Figure EK C ))- In contrast, we 
cannot see extended NW clump in the cycle 0 data which is evident in the SV data (Figure 0Kd)). 
By combining both data (Figures 0K a ))> we can significantly improve source structure by recovering 
both extended and compact features. 


Since we have carried out three epochs of monitoring observations at band 7 in the spectral line 
modes as well as a single-epoch continuum observation, we compared these four observational re¬ 
sults. For the spectral line mode, we notice an apparent difference in the synthesize images between 
epoch to epoch. Because such a difference is not significant for the most compact condensation 
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Fig. 1.— Full uv coverage of ALMA observations, (a) Band 6 data including both compact (SV) 
and extended (cycle 0) configurations, (b) Band 7 data including all of the three epochs in spectral 
line mode and one in continuum mode. The circles indicate the uv length of 100 kA. 


Table 2. Continuum sources at band 6 


ID 

HKKH 

a(J2000) 
+05h35m (s) 

5(J2000) 

-05d22'(") 

Convolved size 
("x ") 

PA 

(deg) 

Integrated 

(mJy) 

Peak 

(mJy beam -1 ) 

1 

14.5373(13) 

27.072(15) 

0.82(3) x 0.61(4) 

127(3) 

123(8) 

74(5) 

2 

14.5988(26) 

28.876(39) 

0.73(6) X 0.48(8) 

110(4) 

123(16) 

106(13) 

3 

14.5902(27) 

29.369(31) 

0.88(6) x 0.52(8) 

97(50) 

160(21) 

104(13) 

4 

14.5073(10) 

30.603(12) 

0.66(2) x 0.47(3) 

124(2) 

284(14) 

272(13) 

5 

14.5804(26) 

31.029(31) 

0.61(6) x 0.57(8) 

97(31) 

118(15) 

102(13) 

6 

14.2775(20) 

30.776(23) 

0.68(6) X 0.46(5) 

170(4) 

53(5) 

51(5) 

7 

14.4853(20) 

32.338(23) 

1.09(5) x 0.51(6) 

121(3) 

267(25) 

141(13) 

8 

14.4162(66) 

33.684(77) 

0.32(17) x 0.20(19) 

74(29) 

13(4) 

58(18) 

9 

14.4360(78) 

34.488(91) 

0.47(18) x 0.35(24) 

99(19) 

27(10) 

49(18) 

10 

14.1939(16) 

34.289(19) 

0.65(4) x 0.44(5) 

100(4) 

69(5) 

71(6) 

11 

14.1112(18) 

36.383(21) 

0.66(4) x 0.48(5) 

134(5) 

57(5) 

52(4) 


Note. — Numbers in parenthesis represent fitting errors determined by the CASA task imfit in unit of 
the last significant digits. 
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Fig. 2.— Continuum emission maps for band 6 data. Synthesized beam size is indicated 
at the bottom-left corner in each panel. The contours start at -5 a level with an interval of 
10a (5, 15, 25, ...). The (0,0) position is the tracking center position of the cycle 0 data, 
RA(J2000)=05h35ml4.1250s, Decl(J2000)=-05d22 , 36.486". The HKKH source ID numbers are 
indicated in panels (a) and (b). (a) Combine both the SV and cycle 0 data. The noise level (la) 
is 9 mJy beam -1 , (b) Same as (a) but excluding those with uv length less than 100 kA. The noise 
level (la) is 5 mJy beam -1 , (c) Cycle 0 data. The noise level (la) is 5 mJy beam -1 . A circle 
indicates a primary beam size of the 12 m antenna, 25 ,/ . (d) SV data. The noise level (la) is 
19 mJy beam -1 . A circle indicates a primary beam size of the 12 m antenna, 25 ,/ . Note that the 
pointing center is A" northeast of that of Cycle 0 data (see panel (c)). 
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Fig. 3.— Continuum emission maps for band 7 data. Synthesized beam size is indicated 
at the bottom-left corner in each panel. The contours start at -5a level with an interval of 
10a (5, 15, 25, ...). The (0,0) position is the tracking center position of the cycle 0 data, 
RA(J2000)=05h35ml4.1250s, Decl(J2000)=-05d22 , 36.486". The HKKH source ID numbers are 
indicated in panels (a) and (b). (a) Combine all of the four epochs including both in spectral line 
mode and continuum mode. The noise level (la) is 9 mJy beam -1 , (b) Same as (a) but excluding 
those with uv length less than 100 kA. The noise level (la) is 5 mJy beam -1 , (c) A single epoch 
data observed in the continuum mode. The noise level (la) is 9 mJy beam -1 . A circle indicates a 
primary beam size of the 12 m antenna, 18 w . (d) Combine three epochs of observations in spectral 
line mode. The noise level (la) is 12 mJy beam -1 . A circle indicates a primary beam size of the 
12 m antenna, 18 /; . 
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identified as Source I, we attribute the reason to the difference in the uv coverage. When com¬ 
pared a image from the single-epoch continuum observation (Figure 01c)) with that of combined 
image of all three epochs observed in the spectral line mode (Figure Eld)), both results are in good 
agreement with each other. Taking into account these effects, we estimate the accuracy of the flux 
measurement to be 20 %. 

We determine positions and intensities of the compact cores by fitting the two-dimensional 
Gaussian on the images made with uv coverages of >100 kA (Figures [21(b) and 01b)) by using the 
CASA task imfit. The flux densities are measured by employing the images after correction of 
primary beam attenuation. Because there could be false detection due to sidelobes or noises, we 
conservatively regard the sources as real detections if their signal-to-noise ratio greater than 5 for 
both bands 6 and 7. As a result, we identify total 11 sources as listed in Tables 0] and 01 In this 
paper, we define the source name as HKKH (an acronym of the authors’ initials) followed by the 
ID number. 

For the band 6 images, the missing flux in the highest resolution images with uv coverage of 
>100 kA (Figure [2](b) ) compared with the image with full-uv coverage (Figure [2](a)) are larger than 
90% for sources HKKH1, HKKH 6 , and HKKH9. In contrast, more than 90% of flux densities at 
band 6 are recovered for sources HKKH2, HKKH4, and HKKH10, implying compact structures. For 
the band 7 images, a fraction of the missing flux in the highest resolution images with uv coverage 
of >100 kA (Figure [3(b)) is less than 50% in comparison with those of full-uv coverage (Figure 
01a)). This difference in the degree of missing flux between bands 6 and 7 would be attributed to 
the different uv coverages and beam sizes. 

For band 7, source HKKH1, HKKH2, and HKKH3 are outside of the primary beamsize of the 
ALMA 12 m antenna while they are within the field of view of band 6 . Although flux calibration is 
not reliable for these sources, we regard them as detection both in bands 6 and 7. Similar to these 
sources, source HKKH5 is detected at the edge of the field of view of band 7 images. In addition, 
the BN object, a famous massive YSO in this region, is outside of the field of view of bands 6 and 
7, as discussed later. 


Discussion 


4.1. Comparison with previous observational results 


There have been number of high-resolution observations of Orion KL with continuum emis¬ 
sions in wide range of wavelengths. We compare some of the highest resolution data reported 
to date. Examples are shown in Figure 01 Positions of most of our detected sources are coinci¬ 
dent with previous inter ferometric maps of the highest resolution millime ter continuum em ission 
(jFriedel fc Weaveii 1201 il l, molecular lines of methyl formate, HCOOCH 3 ilFavre et al.l 1201 ll ) . and 


radio continuum emission (Felli et al 


1993al fbl) as shown in Figure 0] and Table 01 
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Table 3. Continuum sources at band 7 


ID 

HKKH 

q( J2000) 
+05h35m (s) 

<5(J2000) 

-05d22'(") 

Convolved size 
("x ") 

PA 

(deg) 

Integrated 

(mjy) 

Peak 

(mjy beam -1 ) 

1 

14.5372(17) 

26.939(24) 

0.93(5) x 0.66(5) 

104(5) 

505(58) 

156(18) 

2 

14.6141(17) 

28.816(24) 

0.65(5) x 0.57(5) 

33(3) 

591(66) 

299(33) 

3 

14.5844(23) 

29.317(42) 

0.44(10) x 0.33(7) 

131(7) 

161(25) 

212(33) 

4 

14.5113(5) 

30.571(7) 

0.50(2) x 0.44(1) 

126(5) 

1001(34) 

849(29) 

5 

14.5797(21) 

30.852(30) 

0.69(6) x 0.58(6) 

16(12) 

639(91) 

301(43) 

6 

14.2827(09) 

30.801(13) 

0.91(3) x 0.54(3) 

10(2) 

487(29) 

183(11) 

7 

14.4837(18) 

32.238(26) 

0.83(5) x 0.51(5) 

121(4) 

534(64) 

237(29) 

8 

14.4297(16) 

33.582(24) 

0.62(5) x 0.51(5) 

141(8) 

396(44) 

238(26) 

9 

14.4323(29) 

34.440(41) 

0.63(8) x 0.63(9) 

69(10) 

281(55) 

135(26) 

10 

14.1938(10) 

34.116(14) 

0.75(3) x 0.47(3) 

72(2) 

336(22) 

176(11) 

11 

14.1131(16) 

36.327(23) 

0.70(5) x 0.56(5) 

165(7) 

271(29) 

129(14) 


Note. — Numbers in parenthesis represent fitting errors determined by the CASA task imfit in unit of 
the last significant digits. 


Table 4. Identification of continuum sources 


ID 

HKKH 

7 339 GHz 
(mjy beam 1 ) 

J 245 GHz 
(mjy beam -1 ) 

7 90 GHz a 
(mjy beam -1 ) 

Note 

1 

156(18) 

74(5) 

4.89(29) 

C13 

2 

299(33) 

106(13) 

4.39(43) 

C19, MF10? 

3 

212(33) 

104(13) 


C33?, MF10? 

4 

849(29) 

272(13) 

50.09(97) 

C20, MF11?, Source I 

5 

301(43) 

102(13) 

7.08(94) 

C18 

6 

183(11) 

51(5) 

4.18(9) 

C23 IRc7 

7 

237(29) 

141(13) 

6.05(60) 

C21, MF6, SMA1 

8 

238(26) 

58(18) 

5.31(75) 

C22 

9 

135(26) 

49(18) 


C34?, MF2 

10 

176(11) 

71(6) 


MF12, IRc4 

11 

129(14) 

52(4) 


C32, MF1, Source R 


a See Table 1 in lFriedel fe Weaver! 1201 il l. 


Note. — Numbers in parenthesis represent the errors in unit of the last significant dig- 
its. MF and C in colu mn 5 are methyl formate (iFavre et al.l l201lh and 3 mm continuum 
dFriedel W eaver 20 ll]) peaks, respectively. 
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We compare spatial structure of our mi llimeter /submillimete r maps with the mid-infrared 
emission observed with the Subaru telescope ( Okumura et ah| 2011). As can be seen in Figure 21 
the mid-infrared emission shows an anti-correlation with those of millimeter/submillimeter. It is 
easily interpreted that the millimeter/submillimeter emission are emitted predominantly from dust 
cloud with high opacity even at mid-infrared bands. 


We could not detect some of possible counterparts of the millimeter/submillimeter continuum 
sources detected in other tracers (Figured]). There are possible reasons for such differences. First, 
we employ a rather conservative detection criterion that the continuum sources are regarded as 
detection if their peak intensities greater than 5 times the noise level. This is because we eliminate 
contributions from strong sidelobes as possible. As a result, we might miss weak emission in our 
map. Furthermore, we employ the uv lengths longer than 100 kA to search for only compact sources. 
Thus, our continuum images significantly resolves out weak extended emission sources detected in 
the shorter-baselines in previous interferometer observations. 


4.2. Graybody fitting of spectral energy distributions (SED) 


As tabulated in Table [U we identify 11 continuum sources detected both in bands 6 and 
7 images with the uv coverage of greater than 100 kA. It is known that the SED of millime¬ 
ter/submillimeter continuum emission associated with YSOs can be well fitted by the optically 
thin graybody function as a result of thermal dust emission. Thus we fit the observed flux densities 
at 245 GHz, and 339 GHz as listed in Table [4] to the graybody function; 



Aftot 

~W 




( 1 ) 


K„ 



(2) 


where M tot is the total mass, D the distance, k u the total mass opacity, and B u (Td us t) the Planck 
function at the dust temperature of Tdust- The total mass opacity n is expressed as a power law 
function of the freq uency, and k vo = 0.1 cm 2 g _1 at uq = 1.2 THz, or wavelength of 250/xm, is 
adopted (Hi ldebrandlll983l ) assuming the gas to dust mass ratio of 100. If high resolution millimeter 
continuum data at 90 GHz are available (Friedel & We aver 2011), we also employ these results in 
the SED fitting. 

The results of the graybody fitting are summarized in Table [H The power law index of the 
dust opacity, (3 is fitted by using two or three data points as listed in Table [H Because the dust 
temperatures cannot be determined from our continuum observations, we ass ume the temperature 
of 100 K as obtained from the molecular line observations (Favre et ah 2011) and lower value of 
20 K (jEisner et a l. 20081 ) for comparison. The derived mass is anti-correlated with the temperature 
as a function of (exp (hv/kT) — 1), and hence, the lower temperature values would yield an upper 
limit of the mass. The dust opacity index, f3 is close to 1 for most of the sources, except for Source I. 
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The SED of Source I cannot be fitted by a simple optically thin thermal dust emission as discussed 
later. 


The origin of compact continuum sources are attributed to either star-forming or starless cores, 
although detailed nature of each source is still unknown (see discussion about individual sources). If 
the source is associated with embedded YSO, the compact structure of 200-300 AU and high density 
of 10 8 -10 9 cm -3 are indicative of circumstellar disks as an one of the possible origins of millime¬ 
ter/submillimeter sources. The derived masses (0.08Mq-0.27Mq for assumed temperature of 20 K) 
correspond to the smaller end of the masses for disk candidates in another massive star-forming 
region NGC6334I(N) (Hunter et al. 2014). Due to the closer distance of Orion KL than that of 
NGC6334I (N), 1. 3 kpc, derived source sizes are significantly smaller than those in NGC6334I(N) 
(Hunter et al.ll2014l ). On the other hand, the deriv ed mass range is smal ler than the typical circum¬ 


stellar disk masses associated with massive YSOs (Gesaroni et al 


200711 bu t is larger than those of 


low-mass YSOs including members of Orion Nebula Cluster ((Eisner et al. 2008 ). The smaller cir¬ 
cumstellar mass compared with another massive YSO was also recognized for Source I in Orion KL 
by the previous SMA observation (Beuther et al. 2004). If our millimeter/submillimeter sources 
trace the circumstellar disk of massive YSOs, the lower-mass rang e could i mply a relatively later 
evolutionary stages in massive disk formation as proposed in Beuther et al. (2004)- Nevertheless, 
we cannot rule out possibilities of an envelope around YSO, shocked dense gas, or starless dense 
gas, which are significantly resolved with our high-resolution images, or a disk associated with low- 
/intermediate-mass YSO. Key parameters to distinguish these possibilities are gas temperature 
and velocity structure traced by molecular lines in multiple transitions to investigate physical and 
dynamical properties of the continuum sources. 


5. Individual Sources 


5.1. Source I (HKKH4) 


Source I is thought to be a dominant energy source in Orion KL ( Mente n & Reid 1 995 ). 
Number of observational studies have been done for Source I to reveal its nature whereas it can 
be detected only at longer wavelength than submillimeter band due to an extremely high opacity 


(Greenhill et al, 2004; Okumura et al. 


2011 


Sitarski et al 


2013 1. It can be identified as a color 


temperature peak derived from mid-infrared image s (lOkumura et al . 2011|). Position of Source I 
corresponds to the 3 mm continuum source C20 (Friedel fc Weaver! l201ll ). There is a nearb y 
molecular peak MF10 at the western side of the submillimeter emission peak (Favre et al. 2011|). 
Because the position offset between MF10 and^Source I is signi ficant and, in addition, there is no 
emission from organic molecules at Source I (Beuther et al.l 2005), they are unlikely to be physically 
associated. 


By using previous observational results of interferometric continuum observations at centime¬ 
ter, millimeter, and submillimeter wavelengths, we plot a SED of Source I as shown in Figures [5] 
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Fig . 4.— ALMA band 7 continuum maps (contour) superposed on the Subaru mid-infrared image 
(lOkumura et al.ll201lh . Positions of band 7 continuum sources (Table [3]) are in dicated by magent a 
smaller crosses with the H KKH ID numbers. Peak s of methyl formate (MF; iFavre et aL[]2 01lL 


3mm continuum emission (jFriedel &; Weaver] l201lh . and radio continuum emission (IFelli et al 


1993a 


b) are indicated by green squares, blue larger crosses, and yellow triangles, respecti v ely. 


The (0,0) position is that of the 22 GHz HoO maser burst in the Compact Ridge (IHirota et al 


20 11. 2014b). 
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and Table [6] (see Plambeck et ah 2013|, and references therein). At several frequencies, there are 
more than two observational results showing different flux densities. For examp le^flux density of 
our band 6 data, 284 mJy, is consistent with that of CARMA, 310 mJy ( Plambeck et al. 2013|), 
with in the er r ors ^ W%. On the other ha nd, our band 7 data is intermediate between those of 


Beuther et al .1 (j2004f ) and lTang et al .1 (|2010h observed with SMA. The reason for such discrepancies 
could be that insufficient uv coverage affects the flux density measurements due to strong negative 
sidelobes from the bri ght extended emission in the complex structure of the Orion KL region as sug¬ 
gested in Plambeck et al. (2013). Another possibility could be a different degree of missing fluxes 
due to the different uv coverages. To avoid such artifact, we employ the flux density observed 
with the highest resolution observations at each frequency. The data employed in the following 
discussion are indicated by Y in Table [6l 


5.1.1. H free-free emission 


As shown in Figure [5] and Table 0 the SED of Source I can be fitted by a single power law 
function 


log Ff l = p + q log v 


(3) 


with the power law index q of 1.97T0.10 for the frequency range between 8.4 GHz and 690 GHz. The 
power law index is consistent with that expected for the black body radiation, q=2 .0, suggesting 
that the emission could be optically thick even at the submillimeter wavelength. As proposed 
for typical massive YSOs, the emission mechanism of the centimeter to submillimeter continuum 



In order to explain observed SED of Source I from centimeter to submillimeter wavelength, 
we first evaluate physical properties of Source I in the case of the H _ free-free radiation. De- 
taile d derivati on of the H~ free-free radiation is summarized in Appendix following the discussion 
by Reid & Mentenl (1997). Hereafter we assume that Source I is an edge-on disk with uniform 
density an d te m pera ture for simplicity. The apparent size of Source I of 0.20 ,/ x0.03 ,/ is employed 
fPlambeck et ah 2013), corresponding to the disk diameter and thickness of 84 AU and 12.6 AU, 
respectively. With these parameters, a total hydrogen density of 10 12 cm^ 3 corresponds to a column 
density of 1.26xl0 27 cm" 2 (with a diameter of 84 AU) and a total mass of 0.20Mq. 
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As suggested from Figure El Source I seems to be optically thick at 245-339 GHz or even 
at 690 GHz. Thus, the turnover frequency at which optical depth of the H~" free-free radiation 
becomes unity is as high as >300-600 GHz (see Appendix). Using the relationship between the 
total hydrogen density (a sum of molecular and atomic hydrogen) and temperature to explain the 
turnover frequency (Figure [9] and Table [9] in Appendix), we can constrain the range of gas density 
and temperature. 


According to the SMA continuum observation s (lTang et a l. 2010 1. the mass of the main dust 
ridge including Source I is estimated to be 2-12M 0 . Favre et ah (2011) also derive a consistent value 
of 5.8M 0 for their continuum source Ca, in which Source I is located at the western edge of the core. 
The total mass of Source I itself is estimated fr om the velocity structure of the rotating disk traced 


by the SiO masers and submillimeter HA) lines (IKim et al.ll2008l : lMatthews et al. 


201C; 


Hirota et al 


2014al ). >7 Mq, infrared spectroscopy ( Testi et al.1 201C ) , 10AT.), or a momentu m conservation law 


suggested by proper motion measurement (jGoddi et ah 2011; Bally et ah 201 lj), 20 Mq. Thus, the 
circumstellar disk mass of Source I would be much lower than that of the central protostar and its 
host core, <ClOM 0 . With this constraint, the total hydrogen density is lower than 5xl0 13 cm -3 as 
indicated in Figure [9] in Appendix. In this case, the gas temperature would be larger than 1200 K 
for the turnover frequency of >300 GHz (Figurc[9]and Table[9]in Appendix). If the gas temperature 
is about 3000 K or higher, the lower gas density of ~10 n cm -3 is allowed to explain the turnover 
frequency of >300 GHz. 


Using the calculated optical depth, t>h and t v > h 2 for the above temperature (1200-3000 K) 
and density (10 11 -10 14 cm -3 ) ranges, we calculate the flux density F u at each frequency for the 
source size of H; 


F v 


T b 


2u 2 kTs 


dQ 


2iFkT 


(1 -e" 


'^total 


)dn 


(1 - e Ttotal )T 


Total T,H T T.H ‘2 


(4) 

(5) 

( 6 ) 


The results are shown in Figure El Higher temperature or higher density results can explain the 
single power law SED even for the SMA result at 690 GHz. In the case of total hydrogen densities 
of 10 11 cm' 3 and 10 12 cm -3 corresponding to the mass of ~O.O2M 0 -O.2M 0 (Figures EK a) and (b)), 
gas temperatures of higher than 2700 K and 1800 K, respectively, can explain the SED of Source I 
with single power law index of 2.0 up to 339 GHz band. For larger densities of an order of 10 13 cm' 3 
or the mass of a few M 0 (Figure EKc)), the lower temperature of 1500 K can also match the slope 
of the SED of Source I up to 339 GHz. On the other hand, the higher temperature results tend 
to overestimate the flux densities by a factor of ~2. When the emission is optically thick, the 
flux density is proportional to the apparent source size and gas temperature. To account for the 
observed flux density of Source I, the lower gas temperature ~1500 K would be more plausible if 
the source is resolved with the beam size of ALMA. Alternatively, the angular size of Source I, H, 
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would have smaller by a factor of 2 or the source structure would be clumpy with the beam filling 
factor of 50%. Higher resolution observations to resolve the source structure with accurate flux 
calibration will be able to solve this degeneracy of source size and temperature. 


Hirota et al. (l2014al ) reported that the vibrationally excited H 2 O line at 336 GHz is emitted 


from a circumstellar disk. The 336 GHz H 2 O line is thought to be excited thermally with an excita¬ 
tion temperature of >3000 K probably heated via an accretion shock. This excitation temperature 
is slightly higher than that estimated from the present continuum data. According to the velocity 
structure, the 336 GHz H 2 O line map is interpreted as a ring-like structure with a diameter 0.2" 
(84 AU), although the spatial resolution is still insufficient to resolve the structure. It is proposed 
that the ring-like stru ctur e would refl ect t he distribution of the hot molecular (H 2 O) gas probably 
heated via accretion (Hirota et al. 2014a). Alternatively, it is also likely that the 336 GHz H 2 O 
line could trace only the edge of the disk due to the high continuum opacity even at 336 GHz as 
suggested by the SED (Figure [5]). If the optically thick H _ free-free radiation is the main source 
of opacity at 336 GHz, the latter scenario would be more plausible. 


5.1.2. H free-free emission + thermal dust emission 


In Figure 0 our ALMA band 7 data still suggest a marginal excess flux compared with the 
single power law fitting result. The flux density at 690 G Hz also shows a signi ficant excess probably 


origin ated from thermal dust emission as suggested by iBeuther et al .1 ([20061 ) and iPlambeck et al 


(|201 3i). There seems to be systematic deviation between the best-fit model and observed results 
from 43 GHz and 229 GHz, suggesting a change in the emission mechanism at the middle of 
the centimeter to millimeter wavelength. Thus, we next consider the contribution from thermal 
graybody emission from dust. The H _ free-free emission and dust graybody emission contribute to 
the flux densities at lower and higher frequency regions, respectively. In this case, the flux density 
F u at each frequency v is fitted to the sum of power law function due to the H” free-free radiation 
and the graybody radiation from the thermal dust emission; 


Fy = Ff + Ff. 


(7) 


The combined fitting result gives the lower power law index of 1.60T0.24, than that of the 
single-power law fitting, 1.97 (Table [71). The significant amount of the observed flux densities at 
higher frequency (>200 GHz) can be explained by the graybody radiation rather than the H“ free- 
free radiation. In such a case, the turnover frequency could be as low as 200 GHz. The required 
density or temperature under this condition are reduced by a factor of about (2/3) 2 when we change 
the turnover frequency from 300 GHz to 200 GHz (see equation (IA19l) h 

Because of the lack of far-infrared wavelength data corresponding to the flux maximum of the 
graybody function, we cannot constrain the dust temperature of the emitting region. It is possible 
that the dust graybody emission arises from the different volume of gas mainly in a cooler outer layer 
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of Source I or the same volume of gas as the free-free emission at the temperature higher than 
~1000 K. In the former case, the total dust+gas mass contributed by the graybody SED is estimated 
to be 0.082 M 0 assuming the dust to gas mass ratio of 100 and the dust temperature of 100 K (Table 


This value is smaller than that obtained from the SMA observations, 0.2 M 0 (Beuther et al 


20041 1 . In the latter c ase, the dust temperature would be high er than 1000 K as estimated from the 
continuum emission (IReid et al.1 120071 : iPlambeck et al.l 120131 ) and vibrationally excited H 2 O lines 
(Hirota et al. 2014a). The resultant circumstellar dust+gas mass would be reduced by a factor of 
10, ~0.008 M 0 , compared with that of 100 K. Thus, the dust mass would be smaller by a factor 
of 100 (i.e. assumed dust to gas mass ratio), 8xl0 -5 M 0 . This value is much smaller than the gas 
mass estimated from the opacity at this temperature as discussed above (see Figure®. In either 
case, our result may suggest that the dust to gas mass ratio is unusually low in the close vicinity 
of Source I within ~100 AU in diameter because of the dust sublimation under the temperature of 
>1000 K (e.g. Reid et al. 2007 1. 


A power law in dex smaller than 2.0 would reflect a densit y structure of the emitting region 
(e.g. discussion in iBeuther et al.l 120041 : IPlambeck et al.l 1201,81 ). In the present analysis, we do 
not make any correction for above flux densities such as the different source/beam sizes, which 
could more or less induce uncertainties in the derived parameters. To better constrain the emission 
mechanism, higher resolution observations with accurate flux calibration are required to resolve the 
size of Source I. 


5.2. Hot Core (HKKH5) 


Hot Core in Orion KL has been recognized as a prototypical dense and hot molecular gas 
clump which are often found in massive star-forming regions. Hot Core is known to s how wealth 
of molecular line emissions in millimeter and submillimeter wavelengths (e.g. IBeuther et al. 2005 .). 
A submillimeter continuum source SMM3 detected with SMA (Zapata et al- 2011) corresponds to 
this compact millimeter/submillimeter source. This source corresponds to the 3 mm continuum 
emission C18 (iFriedel fe W eaver 20.11) while it cannot be seen in the molecular line map of methyl 
formate (Favre et ah 2011). It is most likely that methyl formate is deficient in the Hot Core 
where nitrogen-bearing organic molecules are dominant. According to the recent molecular line 
observations, Orion Hot Core is most likely h eated externally by the explosive outflow rather than 
by an embedded self-luminous source (Zapata et ah 2011). The mass derived from our ALMA data 
is O.159M 0 at the largest estimate assuming the temperature of 20 K. The dust opacity index 
derived from the graybody SED model is consistent with 1.0. As discussed for Source I, the mass 


of the main dust ridge including Source I and Hot Core is estimated to be 2-12M 0 (ITang et al 


2010) and 5.8M 0 (jFavre et ah 2011). Although the above estimated mass could not resolve more 
compact structures such as Source I and SMA1, our result only recovers a part of the total mass 
of the dense gas. 
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Table 5. graybody fitting results for continuum sources 





20 K 




100 K 



ID 

Diameter 

M 


n(H 2 ) 

A(H 2 ) 

M 


n(H 2 ) 

A(H 2 ) 

HKKH 

(AU) 

(Mq) 

/3 a 

(cm -3 ) 

(cm -2 ) 

(Af 0 ) 

/3 a 

(cm -3 ) 

(cm -2 ) 

1 

313 

0.077(3) 

0.77(2) 

8.6 10 s 

4.0 10 24 

0.009(1) 

0.59(5) 

1.0 10 8 

4.7 10 23 

2 

252 

0.274(35) 

1.30(7) 

5.8 10 9 

2.2 10 25 

0.031(2) 

1.12(4) 

6.6 10 8 

2.5 10 24 

3 

213 

0.082 

0.58 

2.9 10 9 

9.2 10 24 

0.008 

0.27 

2.8 10 8 

9.0 10 23 

4 

215 

0.156(119) 

0.19(40) 

5.4 10 9 

1.7 10 25 

0.018(12) 

-0.00(37) 

6.2 10 8 

2.0 10 24 

5 

257 

0.159(57) 

0.92(19) 

3.2 10 9 

1.2 10 25 

0.018(5) 

0.74(16) 

3.6 10 8 

1.4 10 24 

6 

263 

0.086(55) 

0.90(33) 

1.6 10 9 

6.4 10 24 

0.010(6) 

0.71(30) 

1.9 10 s 

7.4 10 23 

7 

293 

0.175(78) 

0.99(23) 

2.4 10 9 

1.1 10 25 

0.020(10) 

0.80(26) 

2.7 10 8 

1.2 10 24 

8 

158 

0.103(86) 

0.89(43) 

8.9 10 9 

2.1 10 25 

0.012(9) 

0.70(40) 

1.0 10 9 

2.4 10 24 

9 

212 

0.168 

1.51 

6.0 10 9 

1.9 10 25 

0.016 

1.19 

5.7 10 s 

1.8 10 24 

10 

237 

0.145 

1.16 

3.7 10 9 

1.3 10 25 

0.014 

0.87 

3.6 10 s 

1.3 10 24 

11 

249 

0.107 

1.19 

2.4 10 9 

8.8 10 24 

0.010 

0.87 

2.2 10 s 

8.2 10 23 


a Dust opacity index in equation 0- 


Note. — Numbers in parenthesis represent fitting errors (lcr) in unit of the last significant digits. If the 90 GHz data is 
not available, fitting errors cannot be estimated. 



Fig. 5.— Spectral energy distribution (SED) of Source I. A solid blue line indicates the best fit 
single power law model, F v = pu q with the index of q = 1.97 ± 0.10. A solid red line indicates the 
combination of power law and graybody SED. For the power law and graybody terms are shown 
by the red dotted and d ashed lines, respectively. Open squares represent all the data including the 
present ALMA results (Plambeck et ah 2013, and references therein). Filled circles represent the 
data employed in the fitting, as indicated Y in Table [6] 
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Fig. 6.— SED of Source I predicted from H - free-free emission with total hydrogen densities, 
n(H)+2n(H2) of (a) 10 11 cm -3 , (b) 10 12 cm~ 3 , (c) 10 13 cm -3 , and (d) 10 14 cnr 3 . The path length 
is assumed to be 84 AU, which is equal to the major axis of Source I (0.2"). Calculations are done 
with temperatures from 1200 K to 3000 K with a step of 300 K as shown by solid lines. Filled and 
open squares are the same as Figure [5j Note that the SED at the temperature of 1200 K with a 
total hydrogen density of 10 11 cm^ 3 (a) is not plotted because of the lower turnover frequency than 
the plotting range. 
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Table 6. Flux densities of Source I 


Frequency 

[GHz] 

Flux 

[mjy] 

HPBW 

[mas] 

Array 

Epoch 

SED a 

Reference 

8.4 

1.1±0.2 

220 

VLA 

1994.3 

Y 

Menten &; Reid (1995) 

8.4 

1.2±0.1 

262x217 

VLA 

2006.4 


Gomez et al. (2008) 

15 

1.54±0.18 

140x130 

VLA 

1986.3 

Y 

Felli et al. (1993a) 

15 

1.6±0.4 

150 

VLA 

1990 


Felli et al. (1993b) 

22 

5.7±0.9 

109x97 

VLA 

1991.5 

Y 

Forbrich et al. (2008) 

43 

13=1=2 

250 

VLA 

1994.3 


Menten & Reid 119951 

43 

10.8=1=0.6 

1960x1410 

VLA 

1994.9 


Chandler & Wood (1997) 

43 

13 

41x28 

VLA 

2000.9 

Y 

Reid et al. (2007) 

43 

14.5=1=0.7 

170x150 

VLA 

2007.9 


Rodriguez et al. (2009) 

43 

11=1=2 

58x39 

VLA 

2009.0 


Goddi et al. (2011) 

86 

34=1=5 

1000x380 

BIMA 

1995.0 


Plambeck et al. (1995) 

90 

50=1=5 

400x350 

CARMA 

2009.0 

Y 

Friedel & Weaver 120111 

229 

310=1=45 

150x130 

CARMA 

2009.1 

Y 

Plambeck et al. (2013) 

245 

272=1=13 

630x470 

ALMA 

2012 


Present study 

339 

849=1=29 

460x410 

ALMA 

2012 

Y 

Present study 

341 

1400=1=100 

800x700 

SMA 

2009.1 


Tang et al. (2010) 

348 

320=1=48 

780x650 

SMA 

2004.1 


Beuther et al. (2004) 

690 

6700=1=3200 

1400x900 

SMA 

2005.1 

Y 

Beuther et al. (2006) 


Note. — See lPlambeck et al.l J201‘i'l . 
a Y indicates the data employed in the SED fitting. 


Table 7. SED of Source I 


Model 

p 

9 

M(100 K)[M 0 ] 

P 

Single power law 
power law+graybody 

-1.99(19) 

-1.51(30) 

1.97(10) 

1.60(24) 

0.082(49) 

1.34(83) 


Note. — Dust temperature of 100 K is assumed. 

Numbers in parenthesis represent fitting errors (1 <t) in unit of the last significant 
digits. 
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5.3. SMA1 (HKKH7) 


A submillimeter continuum source SMA1 is first identified by Beuthe r et ah (j2004j) hi the 
SMA observations at 348 GHz. This source is identified to be the 3 mm continuum source C21 


(IFriedel fc Weaveiil201ll l and the methyl format e peak MF6 (IFavre et al.ll201lf ). SMAl is proposed 
to be a powering source of the explosive outflow f Beut her fc Ni sscn 20 081), although the origin of th is 


explosive outflow is still under debate fe.g. IZapata et al.ll20D9l : iBallv et al.ll201ll : lGoddi et al.ll2nill ) 
The flux density and peak intensity of our ALMA band 7 data are 534 mJy and 237 nrJy beam -1 
respectively. On the other hand, SMAl was not resolved with the previous SMA observations with 
the flux density and peak intensity of 360 mJy and 360 mJy beam -1 , respectively (IBeuther et al 


2004:). The difference between SMA and ALMA observations could be attributed to the different 
uv coverage. The dust opacity index, f3 is consistent with 1.0. The derived circumstellar mass of 
SMAl is only ~O.18M 0 assuming the dust temperature of 20 K. It is much smaller than the total 
mass of the main dust ridge of this region as discussed above. 


5.4. IRc7 (HKKH6) and IRc4 (HKKH10) 


We clearly see an anti-correlation between the millimeter/submillimeter emission and the mid- 
infrared emission (Figure 0]). However, we detect compact continuum emission sources associated 
with infrared sources IRc4 and IRc7. The nature of these sources whether they are heated internally 
or externally are still unclear. The infrared spectrum of IRc4 observed with the Subaru telescope 
can be fitted to the Planck function with a single temperature of 140 K (Okumura et al. 2011), and 
it is thought to be heated by an external energy source(s) (lOkumura et al. 2011|). On the other 
hand, a color temperature map obtained from longer wavelength observations with SOFIA suggests 
that IRc4 is a self-luminous source (jde Buizer et al.ll2012l L For IRc7, it is proposed that an outflow 


or radiation from Source n would form a fan-like structure in the infrared emission ( Green hill et al 


2004). According to near-infrared polarization observations with Hubble Space Telescope (HST), 


IRc4 is thought to be a reflection nebula illuminated by a nearby source(s) while IRc7 is most 
likely heated by an embedded YSO (jSimpson et al.l 12006). In the 3 mm contin uum map, IRc7 is 
identified as a compact source C23 but IRc4 is only found as a diffuse emission (Friedel & Weaver 


20111 ). We detect compact high density cores with the sizes and H 2 densities of ~200-300 AU and 
>10 8 cm -3 at the positions of IRc4 and IRc7. Circumstellar masses are estimated to be larger than 
O.145M 0 and O.O86M 0 for IRc4 and IRc7, respectively. 


Except for IRc4, IRc7, and BN as discussed below, we cannot detect millimeter/submillimeter 
counterparts at the position of other IRc sources. 




















































- 23 - 


5.5. Compact Ridge (HKKH11) 


We identi fy a compact millim eter/submillimeter emission corr esponding to a submillim eter 
source SMM1 (jZapata et al.l l201ll). 3 mm c ontinuum source C32 (IFriedel &; Weaved l201lh and 
methyl formate peak MF1 (jFavre et al.ll201ll ) in the Com pact Ridge. The m ass of the con t inuum 


source at the Compact Ridge is estimated to be 4AT., (ITang et al.1120101 ) . lEisner et al.1 (120081 1 


identified a compact 1.3 mm continuum source named HC 438 in the Compact Ridge by using 
CARMA at the highest resolution of 0.5 /7 compared with other previous results. They detected 
the 1.3 mm continuum emission with a flux density of 67.8A14.2 mJy. This is consistent with our 
ALMA band 6 data at a similar spatial resolution. The circumstellar mass of HC 438 is derived to 
be 0.085 M 0 with the as sumed d ust t emperature of 20 K (Eisner et al. 2008)- Our result, O.1O5AL 0 
agrees well with that of Eisner et al.1 ( 20081). If the higher temperature of 100 K derived from the 
methyl formate data (Favre et ah 201 f) is employed, the smaller mass of O.OIOATq is obtained. 


As discussed by Favre et ah (120111 ). there are possible counterpart physically associated with 


the millimeter/submillimeter continuum source; a radio continuum source labeled R (jFelli et al 


1993ai). a molecular peak traced by methyl formate lines, MF1 (Favre et al.1 201 lj), and optical 
source Parenago 1822. In addition, an extremely bright 22 GHz HoO maser source appears in the 
Compact Rid ge at the systemic velocity of ~ 6.9 km s _1 and 7.5 km s -1 . It is sometimes called 
supermaser (| Hirota et al. 2011, 2014b, and references therein). The supermaser is located within 
the millimeter/submillimeter continuum source detected with ALMA. 


Favre et al.l (12011) suggest that the supermaser would be related to the shocked molecular gas of 
the Compact Ridge ( Liu et al1l2002 ). It is also suggested that the maser burst could occ ur as a result 
of an interac t ion wit h outflows and a pre-existing YSO in the Compact Ridge (IGarav et al.lll989l ). 
Hirota et al. (2014bj) speculate that the compact continuum source would play an important role to 
amplify the supermaser which can supply plenty of ambient molecular gas. Near the supermaser, a 
cluster of other 22 GHz H 2 O maser features detected with VLA ([Gaume et al. 1 998 ) are distributed 
around the continuumjsource^ It may imply an embedded YSO as a powering source of the H 2 O 
maser cluster (Hirota et al. 2014b). 


We note that an infrared source IRc5 is located at l /7 east of the ALMA continuum peak 
position and is corresponding to another 3 mm continuum peak, C29 ( Friedel Sz Weaver! l201ll ). 
The infra red spect rum of IRc5 can be fitted with the Planck function at the temperature of 
120 K (Okumura et al.l 2011). Similar to IRc4, IRc5 is thought to be a shocked molecular gas 
heated externally (Okumura et al. 2011). Thus, IRc5 is not physically associated with the millime¬ 
ter/submillimeter compact source detected with ALMA. 
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5.6. BN 


Orion BN object is another interesting sou rce regarding the energy source of this region 


(Galvan-Madrid et al 


with Source I (Gomez et al 


2012 ; 


Plambeck et al 


2008: 


Goddi et al 


2013 ) and a possible count erpart of the close encounter 


2011 


Bal ly et al. 2011) or 9 1 C (Ch at ter)ee & Tan 


20121 ). It is a ma ssive YSO associated with a hypercompact Hu region ioniz ed by a zero age main 


seque nce B star ()Plambeck et al.ll2013l ). According to the ALMA SV data (IGalvan-Madrid et al 
2012) with a spatial resolution of 1.32" x 0.62", a compact continuum emission at 230 GHz is de¬ 


tected with a flux density of 126 mJy by using the baselines longer than 100 kA to filter-out extended 
emission. On the other hand, the higher resolution CARMA observatio n detects the 229 GHz con¬ 
tinuum emission of 240 mJy with the resolution of 0.14" (Plambeck et al. 2013|). A continuum 
emission source associated with the BN object is marginally seen in our ALMA bands 6 and 7 
images outside the primary beam size of the ALMA antenna (Figures [2] and [3]). Flux densities at 
bands 6 and 7 are derived to be 55±4 mJy and 33±6mJy, respectively. When we correct for an 
effect of the primary beam attenuation by a factor of 2, the calibrated peak intensities are 110 mJy 
and 66 mJy, r especti vely. Marginal de tect ion in our ALMA image at band 6 does not contradict to 
the results of Galvan-Madrid et al. (120121 ) although our ALMA results are quite uncertain. Thus, 
we do not consider the BN object as an identified source in the present paper. 


5.7. Source n 


Source n is also detected in the SMA observation at 348 GHz l 

Beuther et al. 

2004 

). It is 

proposed to be a YSO associated with a 

bipolar radio jet ( 

Menten & Reid 

1995 

) and a circumstcllar 

disk traced by a mid-infrared emission ( 

Gheenhill et al. 

2004). Our ALMA images at bands 6 and 7 


show no compact emission at the position of Source n. The upper limit (5cr) of their peak intensities 
are 25 mJy beam -1 for both bands 6 and 7. The peak intensity derived from the SMA observation, 
300 mJy beam -1 is much larger than our ALMA result. The non-detection of Source n in our 
ALMA images are possibly due to our lower sensitivities for extended emission components. 


6. Summary 


We have carried out millimeter/submillimeter continuum imaging of the Orion KL region by 
using newly constructed ALMA at a resolution of 0.5", corresponding to a linear scale of 200 AU. 
Compact continuum emission sources are detected and 11 sources are identified at both at band 6 
(245 GHz) and band 7 (339 GHz). They include some of remarkable sources in Orion KL such 
as Source I, Hot Core, SMA1, IRc4, IRc7, and Compact Ridge. Their physical properties such as 
size, mass, H 2 n umber density and column densities are discussed by employing published 3 mm 
continuum data (Friedel & Weaver 2011) to construct SED of dust graybody emission. 
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Among these identified sources, SED of Source I, which is a dominant energy source in 
Orion KL, are extensively studied by using previous obse rvational results from centimeter to sub¬ 
millimeter wavelengths (see references in Plambe ck et al.ll2013 ). The SED model with H~ free-free 
emission is presented following the discussion by Reid & Menten (J.1997) to explain the power-law 
index of the SED, 1.97, consistent with an optically thick emission. We introduce a turnover fre¬ 
quency of the H free-free emission to constrain the gas temperature and total hydrogen density 
of the circumstellar disk of Source I. As a result, total hydrogen density of 10 11 -10 14 cm' 3 are 
required to account for the SED with a single power-law index of 2.0 under the temperature of 
1200-3000 K in the case of turnover frequency of ~300 GHz. When we employ the combination of 
dust graybody emission and power-law SED, the turnover frequency would be as low as 200 GHz, 
which reduces the estimated temperature and/or density by a factor of (2/3) 2 . The fitting result 
yields a smaller power-law index of 1.60, suggesti ng a compact size or clumpy structure of th e emis¬ 
sion region unresolved with the present study (jBeuther et al.1 12004J : iPlambeck et al.1 120131 ) . The 
estimated temperature, density and source size are strongly coupled with each other and hence, 
future higher resolution observations with ALMA will be key issues to solve these degeneracy in 
physical properties. 


We are grateful to S. Okumura for providing a Subaru mid-infrared image. This paper makes 
use of the following ALMA data: ADS/JAO.ALMA^2011.0.00009.SV and 2011.0.00199.S. ALMA 
is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together 
with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. 
The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. We thank the staff at 
ALMA for making observations and reducing the science verification data. T.H. is supported by 
the MEXT/JSPS KAKENHI Grant Numbers 21224002, 24684011, and 25108005, and the ALMA 
Japan Research Grant of NAOJ Chile Observatory, NAOJ-ALMA-0006. M.H. is supported by the 
MEXT/JSPS KAKENHI Grant Numbers 24540242 and 25120007. 

Facilities: ALMA. 


A. H free-free radiation 


The free-free transition of electrons interacting with neutral hydrogen atoms and molecules is 
called H _ free-free emis sion. The formulas for th e absorp tion coefficients o f H~ free-free emission 
were first developed by iDalgarno fc Lanel (119661 ). Then, iReid fc Mentenl 1)19971 ) applied the H - 
free-free emission mechanism to explain radio photospheres of Mira-type variables by simplifying 
the formula of H~ free-free absorption coefficients. Taking into account a similarity in observed 
propertie s between Mira-type variables and Source I, the H~ free-f ree emission can be applied to 
So urce I ( Bcuthe r et al.ll2004l : IReid et al .1120071 : IPlambeck et al.ll2013t ). Here we follow the discussion 
by Reid & Mente n (19971 ) to construct a SED of Source I, assuming the source diameter (path 
length) of 84 AU corresponding to the major axis of Source I of 0.2". 
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The H free-free absorption coefficients at 10 GHz for unit electron pressure per hydrogen 
atom, kiogHz,H; and molecule, /viogHz,h 2 > at the temperature T are approximately expressed as 


kioghz,h = ao,H + ai,nT + a2,uT 2 + a^^T 3 , (Al) 

a 0 , H = +3.376 x 10~ 17 , (A2) 

ai >H = -2.149 x 10" 20 , (A3) 

a 2 ,h = +6.646 x 10" 24 , (A4) 

a 3)H = -7.853 x 10" 28 , (A5) 

and 

«ioghz,h 2 = ao,H 2 + ai,H 2 T + c^h^ 2 + «3,h 2 7" 3 , (A6) 

ao,H 2 = +8.939 x 10” 18 , (A7) 

oi,h 2 = -3.555 x 10 -21 , (A8) 

a 2i H 2 = +1.100 x 10~ 24 , (A9) 

a 3 H 2 = -1.319 x IQ’ 28 , (A10) 


respectively, in unit of cm 4 dvn~ 1 ( R,eid fc Mentenlll997l ). The above equat ions are valid under low 
temperature (<3000 K) and in long wavelength (>3 /mi) as discussed in .Reid &: Menten (1997). 
Because the absorption coefficient is proportional to the inverse square of the frequency, u, then 
the absorption coefficients can be written as 


c+H 


o+h 2 


aio GHz,H 

a 10 GHz,H 


/10 GHz\ 2 

V v GHz ) 
/10 GHz\ 2 

V v GHz ) 


KlO GHz.H+enji/cT 


(10 GHz\ 2 
V v GHz ) 


«io GHz,n 2 n e nji 2 kT 


/10 GHz 
\ v GHz 


2 


(All) 

(A12) 


where k is the Boltzmann constant and n e , ?ih> and n-n 2 are number densities of electron, hydro¬ 
gen atom, and hydrogen molecule, respectively. The optical depths for H~ free-free emission for 
hydrogen atoms and molecules are 


total 

= T v ,n + T^H 2 


(A13) 

Tv, H 

= a U) nL = 

«4o GHz,nn e nuLkT ^ 

10 GHz\ 2 

(A14) 

v GHz ) 




flO GHz\ 2 

(A15) 

Tv, H 2 

= O+H-h = 

«io GHz,H 2 n e nH 2 ++7' 

l i/ GHz ) 


where L represents the path length along a line of sight. 


The optical depths, (1A14I) and (1A15I) . are obtained by solving the thermal equilibrium of 
hydrogen molecules, atoms, and ions at the fixed temperature (Cox 20000; 


nr 


»h 2 


/ lirm^kT 

V h? 


exp - — 


D \ U n (Tf 


kTj Uu 2 (T) ’ 


(A16) 
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n H+ 

nu 


n e 


( 2tt m e kT 
( h? 



2 Um(T) 
Uk(T) 


(A17) 


Here, m e and ttt-h are the mass of electron and hydrogen atom, resp ectively, and mn 2 is a re¬ 
duced mass of the hydrogen molecule, 0.504 atomic mass unit (Cox 2000). D and I are dissociation 
energy of hydrogen molecule, 4.478 eV, and ionization energy of hydrogen atom, 13.598 eV, respec¬ 
tively (jCoxI 120001 ). Partition function of the hydrogen ion and atom are denoted as Un+(T)=1 and 
Uu(T)=2 , respectively, The partition function of the hydrogen molecule is calculated as kT/hcB if 
kT/hcB is much larger thanJL with the rotational constant B in unit of cm -1 . For the hydrogen 
molecule, B is 60.85 cm -1 (Cox 2000) and hence, this approximation is valid in the temperature 
of S>88 K. 


As discussed in Reid & Menten (1997), main sources of electron under the temperature lower 
than ~4000 K are abundant metals such as Na, Mg, Al, Si, K, Ca, and Fe while hydrogen ion is a 
dominant source of electron at the temperature higher than ~4000 K. We evaluate the ionization 
degree of Na, Mg, Al, Si, K, Ca, and Fe, along with H, by solving their Saha equations; 



( 2irm e kT 


y- 


h 2 


exp - — 


I 

kT 


2 Uj(T) 
U 0 (T) ' 


(A18) 


assuming the Solar-svsteni chemical abundances ( 

Cox 

200C). The ionization energy and partition 

function for each species are taken from 

Cox 

(2000) anc 

1 Gray 

to 

o 

o 

Oi 

), respectively. First, the Saha 


equation of each species is solved individually to estimate an ionization degree of each species 
and a total electron density. Then using these initial values, we solve the Saha equations of all 
species iteratively until all of the solutions are converged. In Table [8] we tabulate the calculated 
ionization degree for the temperature range from 1000 K to 3000 K, in which the approximation 
of the H~ free-free absorption coefficients (equations 1 All and IA6I) are valid, with a step of 100 K. 
The results are also shown in Figure (7)(a). Molecular hydrogen is converted to atomic hydrogen at 
the temperature higher than ~2000 K. On the other hand, thermal ionization of hydrogen is not 
significant until the temperature exceeds ~8000 K. 


Using the number densities of electron, atomic and molecular hydrogen, we calculate the 
optical depths of the free-free radiation. Figure 0b) shows the optical depths of the H _ free- 
free radiation as a function of temperature. The dominant opacity source of H~ free-free emission 
is molecular hydrogen under the temperature of 1700 K, and the optical depth due to molecular 
hydrogen shows the maximum at the temperature of 1900 K. At the higher temperature range, 
atomic hydrogen becomes the dominant source of opacity of H~ free-free emission. 


The optical depth is proportional to the inverse square of frequency (Reid fe Ment en 19971 ). 
Even if the H~ free-free radiation is optically thick at lower temperature, the emission becomes 
optically thin at higher frequency band. Thus, we can define the transition frequency from optically 
thick to thin regime as a turnover frequency, u to , as employed in the proton-electron free-free 
radiation. At the turnover frequency, the optical depth becomes unity. From the equations (1 A13D . 
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(| A14f) . and (I A15|) . the turnover frequency can be expressed as 

"to = 10 LkT (KlOGHz,H?LAH + K 10GHz,H 2 n e n H 2 ) GHz. (A19) 

Figures [HJ^a) and [8](b) show the turnover frequency as functions of temperature and total 
hydrogen density, respectively. Although the turnover frequency cannot be expressed by a simple 
formula as a function of temperature, it can be estimated as a simple power law function of the 
total hydrogen density at a given temperature; 

log i/ to = A + B log[n(H) + 2n(H 2 )] (A20) 

The fitting results are summarized in Figure[8jb) and Table[9l Using these results, we constrain 
the gas temperature and the total hydrogen density at the fixed turnover frequency of 200 GHz, 
300 GHz, and 600 GHz as shown in Figure [9] and Table [9l 
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Temparature [K] Temperature [K] 


Fig. 7.— (a) Calculated number density of hydrogen atom, hydrogen molecule, and electron with 
the assumed total hydrogen density of n(H)+2n(H2)=10 12 cm -3 . Electron density is calculated 
from the Saha equations of Na, Mg, Al, Si, K, Ca, Fe, and H (see text), (b) Calculated optical 
depth ofH- free-free emission at 300 GHz. Contributions from H atom and H 2 molecule are plotted 
separately. The total hydrogen density of n(H)+2n(H2)=10 12 cnr 3 and path length of 84 AU are 
assumed. 



Fig. 8.— (a) Turnover frequency as a function of temperature for the fixed total hydrogen density 
from 10 11 cm -3 to 10 15 cm -3 , (b) Turnover frequency as a function of total hydrogen density 
for the fixed temperature from 1200 K to 3000 K with a step of 300 K. Filled and open symbols 
represent the calculated values. Solid lines represent the best fit power-law function as summarized 
in Table [9j The path length is fixed to be 84 AU for both panels. 








Table 8. Fractional abundances of atomic hydrogen, molecular hydrogen and electron 


n k ) 


10 11 cm 3 



10 12 cm 3 



10 13 cm 3 



10 14 cm" 3 


H 

h 2 

e 

H 

h 2 

e 

H 

h 2 

e 

H 

h 2 

e 


1.0(-05) 

5.0(-01) 

2.6(-10) 

3.2(-06) 

5.0(-01) 

5.2(-ll) 

1.0(-06) 

5.0(-01) 

i-3(-ii) 

3.2(-07) 

5.0(-01) 

3.9(-12) 

1100 

1.1 (-04) 

5.0(-01) 

1.4(-09) 

3.4(-05) 

5.0(-01) 

4.2(-10) 

1.1 (-05) 

5.0(-01) 

1.3(-10) 

3.4(-06) 

5.0(-01) 

4.1(-11) 

1200 

8.0(-04) 

5.0(-01) 

8.9(-09) 

2.5 (-04) 

5.0(-01) 

2.9(-09) 

8.0(-05) 

5.0(-01) 

9.5(-10) 

2.5(-05) 

5.0(-01) 

3.0(-10) 

1300 

4.3(-03) 

5.0(-01) 

4.0(-08) 

1.4(-03) 

5.0(-01) 

1.4(-08) 

4.3(-04) 

5.0(-01) 

4.9(-09) 

1.4(-04) 

5.0(-01) 

1.6(-09) 

1400 

1.8(-02) 

4.9(-01) 

l.l(-07) 

5.8(-03) 

5.0(-01) 

5.1(-08) 

1.8(-03) 

5.0(-01) 

1.9(-08) 

5.8(-04) 

5.0(-01) 

6.5(-09) 

1500 

6.2(-02) 

4.7(-01) 

2.0(-07) 

2.0(-02) 

4.9(-01) 

1.1 (-07) 

6.4(-03) 

5.0(-01) 

5.5(-08) 

2.0(-03) 

5.0(-01) 

2.1(-08) 

1600 

1.7(-01) 

4.1(-01) 

4.9(-07) 

5.9(-02) 

4.7(-01) 

2.2(-07) 

1.9(-02) 

4.9(-01) 

1.2(-07) 

6.1(-03) 

5.0(-01) 

5.4(-08) 

1700 

3.9(-01) 

3.0(-01) 

1.1 (-06) 

1.5(-01) 

4.3(-01) 

4.9(-07) 

5.0(-02) 

4.8(-01) 

2.2(-07) 

1.6(-02) 

4.9(-01) 

l.l(-07) 

1800 

6.8(-01) 

1.6(-01) 

1.9(-06) 

3.1(-01) 

3.4(-01) 

1.0(-06) 

l.l(-Ol) 

4.4(-01) 

4.4(-07) 

3.7(-02) 

4.8(-01) 

2.0(-07) 

1900 

8.8(-01) 

5.8(-02) 

2.6(-06) 

5.5(-01) 

2.2(-01) 

1.8(-06) 

2.3(-01) 

3.9(-01) 

8.8(-07) 

7.9(-02) 

4.6(-01) 

3.7(-07) 


9.7(-01) 

1.7(-02) 

3.7(-06) 

7.8(-01) 

l.l(-Ol) 

2.4(-06) 

4.0(-01) 

3.0(-01) 

1.5(-06) 

1.5(-01) 

4.2(-01) 

6.9(-07) 

2100 

9.9(-01) 

5.1(-03) 

5.3(-06) 

9.1(-01) 

4.3(-02) 

3.3(-06) 

6.1(-01) 

1.9(-01) 

2.2(-06) 

2.7(-01) 

3.7(-01) 

1.2(-06) | 

2200 

1.0(+00) 

1.6(-03) 

6.7(-06) 

9.7(-01) 

1.5(-02) 

4.6(-06) 

7.9(-01) 

1.0(-01) 

2.9(-06) 

4.2(-01) 

2.9(-01) 

1.8(-06) co 

2300 

1.0(+00) 

5.7(-04) 

8.0(-06) 

9.9(-01) 

5.6(-03) 

6.0(-06) 

9.1(-01) 

4.7(-02) 

3.9(-06) 

5.9(-01) 

2.0(-01) 

2.4(-06) ° 

2400 

1.0(+00) 

2.2(-04) 

l.l(-05) 

1.0(+00) 

2.2(-03) 

7.3(-06) 

9.6(-01) 

2.0(-02) 

5.1(-06) 

7.5(-01) 

1.2(-01) 

3.2(-06) 

2500 

1.0(+00) 

9.0(-05) 

1.6(-05) 

1.0(+00) 

9.0(-04) 

9.0(-06) 

9.8(-01) 

8.7(-03) 

6.3(-06) 

8.6(-01) 

6.8(-02) 

4.1(-06) 

2600 

1.0(+00) 

4.0(-05) 

2.4(-05) 

1.0(+00) 

4.0(-04) 

1.2(-05) 

9.9(-01) 

3.9(-03) 

7.7(-06) 

9.3(-01) 

3.5(-02) 

5.2(-06) 

2700 

1.0(+00) 

1.9(-05) 

3.4(-05) 

1.0(+00) 

1.9 (-04) 

1.8(-05) 

1.0(+00) 

1.9(-03) 

9.7(-06) 

9.7(-01) 

1.7(-02) 

6.3(-06) 

2800 

1.0(+00) 

9.2(-06) 

4.3(-05) 

1.0(+00) 

9.2(-05) 

2.6(-05) 

1.0(+00) 

9.2(-04) 

1.3(-05) 

9.8(-01) 

8.9(-03) 

7.7(-06) 

2900 

1.0(+00) 

4.8(-06) 

5.2(-05) 

1.0(+00) 

4.8(-05) 

3.4(-05) 

1.0(+00) 

4.8(-04) 

1.8(-05) 

9.9(-01) 

4.7(-03) 

9.7(-06) 


1.0(+00) 

2.6(-06) 

6.3(-05) 

1.0(+00) 

2.6(-05) 

4.3(-05) 

1.0(+00) 

2.6(-04) 

2.5(-05) 

9.9(-01) 

2.6(-03) 

1.3(-05) 


Note. — Abundances are calculated by n(X)/(n(H)+2n(H2)). Values indicated by a(b) mean a x 10 b . 
The abundance of hydrogen ion is negligible under the condition of low temperature (<4000 K). 
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Table 9. Temperature and total hydrogen density for given turnover frequencies. 


T(K) 

Power law fit a 

n(H)+2n(H2) (cm 

~ 3 ) 

a 

b 

200 GHz 

300 GHz 

600 GHz 

1000 

-8.84 

0.750 

7.1(14) 

1.2(15) 

3.1(15) 

1100 

-8.32 

0.751 

1.4(14) 

2.4(14) 

6.0(14) 

1200 

-7.91 

0.753 

3.6(13) 

6.2(13) 

1.5(14) 

1300 

-7.64 

0.761 

1.2(13) 

2.0(13) 

5.0(13) 

1400 

-7.57 

0.778 

5.0(12) 

8.3(12) 

2.0(13) 

1500 

-7.61 

0.798 

2.6(12) 

4.3(12) 

1.0(13) 

1600 

-7.41 

0.800 

1.4(12) 

2.3(12) 

5.5(12) 

1700 

-7.03 

0.787 

7.2(11) 

1.2(12) 

2.9(12) 

1800 

-6.70 

0.777 

3.9(11) 

6.6(11) 

1.6(12) 

1900 

-6.55 

0.777 

2.4(11) 

4.1(11) 

1.0(12) 

2000 

-6.50 

0.785 

1.7(11) 

2.8(11) 

6.8(11) 

2100 

-6.55 

0.797 

1.2(11) 

2.1(11) 

5.0(11) 

2200 

-6.67 

0.815 

1.0(11) 

1.7(11) 

4.0(H) 

2300 

-6.84 

0.834 

9.0(10) 

1.5(11) 

3.4(11) 

2400 

-6.95 

0.849 

7.9(10) 

1.3(11) 

2.9(11) 

2500 

-6.97 

0.856 

6.6(10) 

1.1(11) 

2.4(11) 

2600 

-6.90 

0.858 

5.4(10) 

8.6(10) 

1.9(11) 

2700 

-6.80 

0.856 

4.3(10) 

6.9(10) 

1.6(11) 

2800 

-6.70 

0.854 

3.5(10) 

5.6(10) 

1.3(11) 

2900 

-6.62 

0.853 

2.9(10) 

4.6(10) 

1.0(11) 

3000 

-6.59 

0.856 

2.5(10) 

4.0(10) 

8.9(10) 


Note. — Total hydrogen densities, n(H)+2n(H2), are calculated 
with fixed turnover frequencies and temperatures. The path length 
is fixed to be 84 AU. 


See equation d A20b . Values indicated by a(b ) mean a x 10 b . 
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Fig. 9.— Relationship between temperature and total hydrogen density, n(H)+2n(H2) for the 
fixed turnover frequency. Filled square, open circle, and filled diamond represents the turnover 
frequency of 200 GHz, 300 GHz, and 600 GHz, respectively. The path length is fixed to be 84 AU. 




